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Abstract 

In this paper, we study synchronization in the delayed discrete-time complex networks. Several criterions 
of synchronization stability for such networks are established. And illustrative examples are presented. The 
numerical simulations coincide with the theoretical analysis. 
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1. Introduction 

Time delays commonly exist in the real networks, some of them are trivial so can be ignorant 
whilst some of can not be ignored, such as in the long-distance communication and traffic congestions, 
etc. So such networks with retard time attracts much attention. Recently, Li, et al, studied the 
following continuous network model with time delays [1], 

N N 

Xi = f(xj) + e}^CjjA(xji(t - n),Xj2(t - r 2 ), • • • ,x jn (t-T n )) T = f(x i ) + ey j c ij A ■ Xj(t - t) , (1) 
i=i 3=1 

where / : 5R n — > 5ft 71 is a continuously differentiable function, X4 — (xn, X42, ■ ■ ■ ,Xi n ) T € K n are 
the state variables of node i, e > represents the coupling strength, A — (a^ ) nx „ € 5ft™ xn indicates 
inner-coupling between the elements of the node itself, while C = (cjj)jVxiV denotes the outer-coupling 
between the nodes of the whole network (it is often assumed that there is at most one connection 
between node i and another node j, and that there are no isolated clusters, i.e., C is an irreducible 
matrix). The entries are defined as follows: if there is a connection between node i and node 
j (j =/= i), then we set cy = 1; otherwise cy = (j ^ i), and the diagonal elements of C are defined by 
<Hi = _ Ylj=i,fr% i = 1, 2, • • • , AT, Tj, i = 1, • • • , n, are the time delays. 

In this paper, we study the discrete version of network (1), which is described as, 



Xi(k + l) = f(Xi(k))+s J2j=l c ijA(Xji(k-Ti),X j2 (k-T2),--- ,x jn {k - T n )) T 



fix^k)) +£^j = i CijA ■ xj(k - t) 



(2) 



in which /, e, Xi (i — 1,2,- ■• , N), C and A have the same meanings as those in (1). The only 
difference is that in Eq. (2), Xi and Tj are defined in the positive integer set Z + . 

In model (2), the inner-coupling is linear since A is a matrix, a natural generalization is that the 
inner-coupling can be nonlinear [2]. Such a discrete-time network with retard time reads as, 



Xi(k + 1) = f(xi(k)) +£ J2f=i CijAg(xji(k - T 1 ),Xj 2 (k - r 2 ), • • • , x jn (k - r„))' i 



f{xi(k)) + eE^Li CijAg(xj(k - t)) . 



(3) 
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where /, e, Xi, t* (i = 1, 2, • • • , N) and C have the same meanings as those in (2), g : 5R™ — > 5R™ is a 
continuously differentiable nonlinear function. 

The outer-coupling configuration C in networks (2) and (3) has following properties [3,4]. 

Lemma 1: Suppose that C = (Cy)jvxJV is a real symmetric and irreducible matrix, where > 

(« ^ j), c l4 = - EjLi,j^» C U' thcn 

(1) is an eigenvalue of C with multiplicity 1, associated with eigenvector (1, 1, • • • ,1) T ; 

(2) all the other eigenvalues of C are less than 0; 

(3) there exists a unitary matrix, $ = (</>i, fa, ■ ■ ■ , 4>n) such that 

C T <j) k = \ k k , k = l,2,---,N, 

where = Ai > A 2 > A3 > • • • > A^v are the eigenvalues of C. 

In what follows, the definition of synchronization for networks (2) and (3) is introduced below [2]. 
Definition 1: Let A be an attractor of the discrete-time dynamical system 

s(k + l) = f(s(k)). 

We say that networks (2) and (3) are (asymptotically) synchronized to A, if for k — ► +00, 

x 'i ^ ^4. - t — - 1 - ■ . -/V". 

In the rest of this paper, the criterions of synchronization stability for networks (2) and (3) are 
established in Section 2. And the numerical examples are presented in Section 3. 

2. Synchronization theorems 

By utilizing Lemma 1, one can derive the following theorem. 

Theorem 1: Consider the delayed discrete-time network (2), let = Ai > A 2 > • • • > Ajy be the 
eigenvalues of the coupling configuration matrix C. If the following N — 1 systems of n-dimensional 
linear delayed equations are asymptotically stable about their zero solutions, 

r](k + 1) = Df(s(k)) V (k) + eXiA ■ n(k - r) , i = 2,---,N, (4) 

where Df(s(k)) £ 5R" X " is the Jacobian of f(x(k)) at s(k), r,{k) G 3? n , n{k - t) = (jyi(fc-n), • ■ ■ ,n n (k- 
t„)) t e SR", s(k) is the orbit of attractor A of equation s(k + 1) = f(s(k)), then network (2) is syn- 
chronized to the attractor A. 

Proof: Linearizing (2) at s(k) yields, 

N 

ei (k + 1) = Df{s(k) ei {k)) +eJ2 c io A ■ ej(fc-r), 1 < i < N, 

where ej(fc) denotes the deviation from the state s(k), i.e., ej(fc) = Xi{k) — s(k), ej(k — r) = (eji(k — 
n), ■ • ■ , ej n (k — T n )) T e 3t n . This linearized system can be rewritten as 

e;(fc + 1) = Df(s(k))ei(k) + sA ■ (ei(k - r) ,e 2 (k - r) , • • • ,ejv(fc - r) )(ca, • ■ ■ ,c iN ) T . 

Let e(fc) = (ei(/c), e2(fc), • • • , ejv(fc)) S 3? TlxJV , the above equation can be expressed in a compact form, 

e(k + 1) = Df(s(k))e(k) + eA ■ ejk - t) C t , 
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where e(k - r) = (ei(fc - r) , e 2 (fc - r) , ■ • ■ ,ejv(fc-r)) e 3?" xjv . 

From Lemma 1, there exists a nonsingular matrix such that C T $ = $r, T = diag(Ai, • • • , Xn)- 
If one sets e(fc)<& = v(k) = («i(fc), w 2 (fc), • • • ,vjy(k)) <E $l nxN , then the above equation can be trans- 
formed into the following matrix equation 



v(k + 1) = Df(s(k))v(k) + eAv(k - r) T, 
which immediately follows that, 



v t {k + 1) = Df(s(k))vi(k) + eXiA ■ Vi {k - r) , i = 1, • • • , N, 



where - r) = (^i(fc - n), • • • ,^„(fc - t„)) t G 5ft™. 

Note that Xi = corresponds to the linearized system of the isolate equation s(k + 1) = f(s(k)). 
If the following TV — 1 pieces of n-dimensional linear time-varying delayed equations 



v i (k + l) = Df(s(k))v i {k) + sX i A-v i (k-T), % = %■■■ ,N, 

are asymptotically stable around their zero solutions, then e(k) will tend to zero, which shows that 
the conclusion holds. This completes the proof. 

Theorem 2: Assume that all eigenvalues of the matrix C in (2) are listed in order, = Ai > A 2 > 
• • • > X N . If there exists a positive-definite matrix P > such that 



Df(s(k)) T PDf(s(k)) -P + I eX l D.f{s{k)) T PA 
eXiA T PDf{s(k)) e 2 X 2 l A T PA - I 



< 0, i = 2,3,- 



,N, 



(5) 



where / is the identity matrix, s(k) is the orbit of attractor A of the equation s(k + 1) = f(s(k)), 
then network (2) is synchronized to A for any fixed delay € Z + (i — 1, 2, • • • , n). 
Proof: Select a Lyapunov functional as 

n k—1 

v( v (k))=v(k) T Pv(k) + Y, E ^(v) T m(*), 

i—1 cr—k — Ti 

in which r/(k) = (r]i(k), r]2(k), ■ ■ ■ ,rj n (k)) T . Then, along the solution of the ith (i = 2,3, ••• ,N) 
equation in system (4), one gets 



AV(r](k)) 



V( V (k + l))-V( V (k)) 



[Df{s(k))ri(k) + eX t A ■ V (k - r)] T P[Df(s{k))r,(k) + eX t A ■ V (k - r) 



-r](k) T Pr](k) + r^{k) T ri{k) - -q{k - r) • r](k - r) 

Df(s(k)) T PDf(s(k)) -P + I eX t Df( s {k)) T PA 
eX t A T PDf(s{k)) e 2 X 2 A T PA - 1 



rj(k) 


T r 


r){k — t) 





rj(k) 
r](k — t) 



, N, which 



By using linear matrix inequality (LMI) (5), one has AV(r](k)) < for all i = 2,3. 
implies that the zero solutions of systems (4) are asymptotically stable. From Theorem 1, network 
(2) is synchronized to A. The proof is finished. 



In the following, one can similarly obtain theorems 3 and 4 for network (3). 
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Theorem 3: Consider the delayed discrete-time network (3), let = Ai > A 2 > • • • > Ajv be the 
eigenvalues of the coupling configuration matrix C. If the following N — 1 systems of n-dimensional 
linear delayed equations are asymptotically stable about their zero solutions, 

rj(k + 1) = Df(s{k))t](k) + e\iADg(s(k - r)) • V (k - r) , i = 2, • • • , N, 

where D f(s(k)), Dg(s(k — r)) 6 3i nxn are the Jacobians of f(x(k)),g(x(k)) at s(k) and s(k — r), 
7](k) G 3£", s(fc - r) - ( Sl (fc - n), ■ ■ ■ , s„(fc - r„)) T e ry(fc - 7) = {m{k - n), • • • , r?„(fc - r„)) T e 
3?™, s(/c) is the orbit of attractor A of the equation s(k+ 1) = f(s(k)), then network (3) is synchronized 
to A 

Theorem 4: Assume that all eigenvalues of the matrix C in (3) are listed in order, = Ai > A 2 > 
■ • ■ > Ajv- If there exists a positive-definite matrix P > such that 

Df{ S {k)) T PDf(s(k))-P + I eXiDf( S (k)) T PADg( S (k-T)) 

< 0, i = 2, • • • ,7V, 

e\Dg(s(k - T)) T A T PDf{s{k)) e 2 X 2 l Dg{s{k - t)) t A t PADg(s(k — t)) — I 

where / is the identity matrix, s(k) is the orbit of attractor A of equation s(k + 1) = f(s(k)), then 
network (3) is synchronized to A for any fixed delay n G Z + (i = 1, 2, ■ • • , n). 

In [1], by using "matrix measure" [5,6], Li, et ai, discussed the synchronization of network (1) and 
the following network (6), 

N N 
Xi =f(x i )+e^2c ij A(x j - L (t-T j - L ),X j 2(t-T j2 ),--- ,X Jn (t-T Jn )) T = f{x i )+£ 1 ^CijA-X j {t - Tj) , (6) 

where /, s, Xi (i = 1, 2, • • • , N), C and A have the same meanings as those in (1), the sole difference 
is that in (6) a different node j has a different time-delay vector, (tji , Tj2, • ■ • , Tj n ). 

We can also apply "LMIs" presented here to establishing synchronization theorems for delayed 
continuous-time networks [1] and [6]. In the following, the definition of synchronization is given. 

Definition 2 [1,2]: Assume that B is an attractor of the continuous dynamical system 

* = /(*)■ 

We say that networks (1) and (6) are (asymptotically) synchronized to B, if for t — ► +oo, 

Xj — ► B, i = 1, ■ • • , N. 



Here, we list the synchronization theorems of networks (1) and (6) just for reference later on. For 
details of the proofs, see [7] . 

Theorem 5: Consider the delayed continuous network (1), let — X\ > A 2 > • • • > Xjy be the 
eigenvalues of the coupling configuration matrix C. If the following N — 1 systems of n-dimensional 
linear time- varying delayed differential equations are asymptotically stable about their zero solutions, 

fj(t) = Df(s(t))n(t) + eX t A ■ t](t - r) , i = 2, • • • , JV, 
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where Df(s(t)) G 3?" x ™ is the Jacobian of f(x(t)) at s(t), rj(t) G K n , r?(t - r) = (771 (i-n), • • • ,r) n (t- 
T n )) T G 5ft", s(i) is an orbit of attractor /3 of equation s = f(s), then network (1) is synchronized to 
B. 

Theorem 6: Assume that all eigenvalues of the matrix C in (1) are listed in order, = \i > X 2 > 
■ ■ ■ > Ajv- If there exists a positive-definite matrix P > such that 



PDf(s(t)) + Df(s(t)) T P + I e\PA 
e\,A T P -I 



< 0, i = 2,3, - - ■ ,N, 



where I is the identity matrix, s(t) is an orbit of the attractor B of the equation s = f(s), then 
network (1) is synchronized to B for any fixed delay > (k = 1, 2, ■ • ■ ,n). 

Theorem 7: Consider the delayed dynamical network (6), let = Ai > A 2 > • • • > \n be the 
eigenvalues of the coupling configuration matrix C. If the following N — 1 systems of n-dimensional 
linear time- varying delayed differential equations are asymptotically stable about their zero solutions, 



fi{t) = Df(s{t))r){t) + eXiA ■ rj{t - n) , i = 2, ■ ■ ■ , N, 



where Df(s(t)) G 5R" xn is the Jacobian of f(x(t)) at s(rj), r?(r;) G K", 77 (i - t,) = (»7i(t— m), • • ■ ,7y„(t- 
7"in)) T G 3?", s(i) is the stable equilibrium B of equation s = f(s), then network (6) is synchronized 
to Bo- 
Theorem 8: Assume that all eigenvalues of the matrix C in (6) are listed in order, — Ai > A2 > 
• • • > Ajv- If there exists a positive-definite matrix P > such that 



PDf(s(t)) + Df{s(t)) T P + I eXiPA 
e\,A T P -I 



< 0, i = 2,3, - - ■ ,N, 



where / is the identity matrix, s(t) is the equilibrium Bq of equation s = f(s), then network (6) is 
synchronized to Bq for any fixed delay tu > (k, I = 1, 2, • • • , n). 

Remark 1: Obviously, Bo C B. B also contains other stable limit sets, e.g., stable (quasi-)periodic 
orbit, strange attractor. Theorems 7 and 8 hold only for Bo- For the other cases, the studies are not 
easy but need further consideration [8] . 



3. Illustrative examples 

In this section, we consider a five-node network, in which each node is a simple 2-dimensional 
Henon map [9,10], described by 



fi(xi,x 2 ) 




1 + x 2 — ax\ 


f2(xi,X 2 ) 




bx 1 



(6) 



If a = 0.5, b = 0.3, (6) has a period solution [9]. 
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We use the coupled configuration matrix (Cy)jvxiv [11], which is 



C 



-2 


1 








I— 1 


I— 1 


-3 


1 


1 








1 


-2 


1 








1 


1 


-3 


1 


I— 1 








1 


-2 



whose eigenvalues are Ai = 0, A 2 = —1.382, A 3 = —2.382, A 4 = —3.168, A 4 = —4.168, and use the 
inner-coupling matrix A = diag(l, 1). 

At first, the time-delay vector (n, t 2 ) = (1,2) is considered. 

Here, set the coupling strength e = 0.015. By using the MATLAB LMI Toolbox, one can obtain 
that there exists a common positive-definite matrix, 



5.3626 
9.1633 



such that the condition of Theorem 2 is satisfied. From Theorem 2, we know that this delayed network 
is synchronized to the stable periodic state of the isolated Henon map. The following numerical 
simulations are also in line with the above theoretical analysis. 

In Fig 1, we plot the curves of the synchronization errors between the states of node i and the 1st 
node, that is, = Xij(k) — x\j(k) for i = 2, • • • , 5,j = 1, 2. 





a) 



40 60 

k 



-0.1 W:— 



b) 



Fig 1. Synchronization errors for the coupled network with time-delay vector (n, T2) = (1, 2). (a) j = 1 (b) j = 2. 

Secondly, we choose the nonlinear inner-coupling function as 



(7) 



By using the MATLAB LMI Toolbox again, we can find there exists a positive-definite matrix 

5.3649 
9.1714 



9i{xi,x 2 ) 






92{xi,x 2 ) 




sin x 2 
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such that the condition of Theorem 4 is satisfied. So the five-node Henon network with time delay 
(t T 2) = (1, 2) is synchronized to the stable periodic state of the isolated Henon map. The numerical 
simulation below shows this point of view. 

In Fig 2, we plot the curves of the synchronization errors between the states of node i and node 
1, i.e., dj = Xij(k) - xij(k) for i = 2, • • • , 5, j = 1,2. 



a) 




k b) 

Fig 2. Synchronization errors for the five-node Henon network, (a) j = 1 (b) j = 2. 
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